Explosive site percolation with a product rule 
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We study the site percolation under Achlioptas process (AP) with a product rule in a 2 — 
dimensional (2D) square lattice. From the measurement of the cluster size distribution, Ps, we 
find that Ps has a very robust power-law regime followed by a stable hump near the transition 
threshold. Based on the careful analysis on the Ps distribution, we show that the transition should 
be discontinuous. The existence of the hysteresis loop in order parameter also verifies that the tran- 
sition is discontinuous in 2D. Moreover we also show that the transition nature from the product 
rule is not the same as that from a sum rule in 2D. 
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The percolation transition describing the emergence 
of large-scale connectivity in lattice systems or complex 
networks has been extensively studied in statistical me- 
chanics and related fields due to its possible applications 
to various phenomena such as sol-gel transition and poly- 
merization, resistor networks, and epidemic spreading [Ij. 
When the occupation probability of node (site) is lower 
than certain threshold pc, all the clusters are microscopic. 
As the occupation probability increases, the macroscop- 
ically connected cluster emerges. Such transition in the 
ordinary percolation is a continuous transition [l[ . 

On the other hand, there have been several attempts 
to find a percolation model which undergoes a discon- 
tinuous transition. The discontinuous percolation tran- 
sition can be found in the modelling of magnetic sys- 
tems with significant competition between exchange and 
crystal-field interactions 0, 0] ■ The similar phenomena 
has been found in financial systems Ql, in which two 
equally probable phases exit. Other examples of the dis- 
continuous transition in percolation are the formation of 
infinite cluster under a central-force |5|1 and the cascade 
of failure in interdependent networks [6| . 

Recently, Achlioptas et al. 0| suggested a simple pro- 
cess in which the growth of large clusters is systematically 
suppressed and the process is usually called as Achlioptas 
process (AP). Based on the analysis of transition inter- 
val it was argued that the percolation transition under 
AP is explosive and discontinuous. Several variant of 
models have been investigated to understand the general 
properties and conditions which cause such non-trivial 
discontinuous transition 0-ll|. Some examples of such 



non-trivial transition has been found in nano-tube based 



system ll| , protein homology network 12| , and commu- 
nity formation [13 1. 

However, more recent studies on the percolation tran- 
sition under AP reveals several evidences which strongly 
suggest that the transition can be continuous. For exam- 
ple, da Costa et al. [iJl argued that the transition in the 
complete graph (CP) is continuous, even though the or- 
der parameter exponent is very small (/3 ~ 0.056). From 
the measurement of the cluster size distribution Lee et al. 



Iq also argued that the transition in CP is continuous. 
Grassberger et al. [3] also argued that the transition, 
even in the low-dimensional systems, can be continuous 
based on a measurement of the order parameter distri- 
bution. 

Since most of the studies on the criticality of AP pro- 
cess are restricted to the infinite dimensional systems, it 
is still not clear whether AP also produces a continuous 
transition in lower dimensional systems or not. For exam- 
ple, in the bond percolation under AP in a 2-dimensional 
(2D) square lattice, the product rule was argued to pro- 
duces a discontinuous transition based on a finite-size 
scaling 0, In contrast Grassberger et al. [3l was 
argued that transition of the 2D AP bond percolation is 
still continuous. We therefore cannot exclude the possi- 
bility that the transition nature of the AP in the mean- 
field limit can be different from that in lower-dimensional 
systems, like the Potts model jlSl . Moreover, based on 
the measurement of hysteresis [l9[ , a sum rule for the 2D 
site percolation possibly makes the transition continuous 
in the thermodynamic limit. This indicates that under 
the AP-like processes the bond percolation and site per- 
colation may have different transition natures in the 2D 
lattice. In the ordinary percolation, bond and site perco- 
lations are known to belong to the same universality class 
[l|. In contrast, the results in Refs. [s, 17, 1^ show the 
possibility that under AP the bond percolation with a 
product rule and the site percolation with a sum rule do 
not belong to the same universality class. Therefore, it 
is theoretically important and interesting to investigate 
whether in a low-dimensional system the product rule 
and the sum rule belong to the same universality class or 
not. In order to achieve this purpose, we investigate the 
site percolation under AP with a product rule and show 
that AP with the product rule produces a clear discon- 
tinuous transition in a 2D lattice. For this we carefully 
analyze the cluster size distribution and hysteresis. 

AP in 2D site percolation is defined as follows: (I) 
We select two sites a and /3 at random. (II) Let 
{sai,Sa2,--- ,Sa„} {{s Pi , s , ■ ■ ■ ,S/3„}) be the sizes of 
clusters which form into a new big cluster with the size 



2 



ELi + 1 (Er=i S/9fc + 1) by occupying the site a 
Here the cluster size is defined by the number of sites in 
the cluster. Then calculate the products 



n 



and 



TTfi 



(1) 



This rule is generally called product rule (PR). (HI) If 
< TT^ (tTq > TT^) then site /3 (a) remains to be va- 
cant. The processes (II) and (III) prefer the connection 
between small clusters, which causes the cluster repulsion 
or suppress the growth of large cluster. If the product in 
Eq. (IlJ is replaced by summation, then the rule is called 
as a sum rule. Recent study for site percolation with sum 
rule shows that the transition is continuous when the lin- 
ear size of the lattice, L, goes to infinity [3]. Since we 
use a 2D square lattice n(rn) in Eq. ([T]) is at most 4. 
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FIG. 1. (color online) (a) Plot of Ps for p = 0.75 or p < pc- 
The data for L = 1024 are denoted by black squares and 
those for L = 2048 are denoted by red circles, (b) Plot of Pa 
for p = 0.8 or p > pc- The solid line represents the relation 
Ps ^ s^^ with 5 ~ 0.90(2). The data points marked by 
"x" represents the Pa for macroscopically large clusters or 
the largest cluster, (c) The same plot for p = 0.77 or near 
the transition threshold, (d) Plot of Ps for p = 0.765 and 
p = 0.770 on the lattice with L = 2048. 

To understand the percolation transition physically to 
a deeper level, the properties of the cluster size distri- 
bution should be the first one to understand [J, [l^ E^l ■ 
The cluster size distribution Ps{p) at a site occupation 
probability p is normally defined by the probability that 
a randomly selected site belongs to a cluster which has 
s sites (s-cluster) . For ordinary percolation, it is well 
known that Ps{Pc) satisfies a power-law relation 



(2) 



with d ~ 1.05 at the percolation transition probability pc- 
Since, as we shall show, the percolation properties under 
AP depend on the history how the clusters are grown, we 
measure Ps{p) by filling sites from the vacant lattice or 



increasing p. In Fig. [l]Ps(p)'s for 2D site percolation un- 
der AP with product rule (2DSAP) are displayed. Ps in 
Fig. [T] shows an anomalously unique behavior compared 
to that of the ordinary percolation (OP) [l| and that 
of the AP percolation on the complete graph (APCG), 
which was argued to undergo the continuous phase tran- 
sition |15 |. 

When p < Pc, Ps for 2DSAP has a hump in the tail 
as p approaches to Pc (l^. In this regime, Pg does not 
depend on L or 7V(= L x L) as shown in Fig. [Tfa). In 
OP, Ps normally decays exponentially as s gets larger 
in this regime. For the detailed comparison to those of 
OP and APCG, let's call s at which the hump is max- 
imal sh- In OP we cannot identify sh- In APCG sh 
and P{sh) was argued to satisfy the scaling behavior, 
iV^ and P{Sh) 



15 



s// ~ iV" and P[f)H) — N~y with a; > and y > Q 
Therefore in APCG the hump has the negligible contri- 
bution and Ps satisfies the same scaling form as Eq. ([2]) in 
the limit N — )■ oo. This Ps behavior in APCG is believed 
to be one of the signals for the continuous transition as 
in OP. In contrast sh and P{sh) of 2DSAP do not de- 
pend on L or Sh — const and P{sh) — const as L{N) 
gets larger. We have numerically checked this behav- 
ior for N = 216,218,220,222. This behavior for 2DSAP 
means that there should be many considerably large sta- 
ble microscopic s-clusters with s ~ sh before transition, 
which indicates the unstable or sudden appearance of the 
macroscopic cluster by connecting these clusters when p 
increases. 

Even when p > pc, Ps for 2DSAP has a unique behav- 
ior as shown in Fig. [Ijb). Except Ps for the macroscop- 
ically large cluster, Ps for microscopically finite clusters 
for p > Pc still satisfies the same power-law Pg = As^^ 
with the same exponent S as Ps{Pc) or S = 0.90(3), which 
we will explain with the data in Fig. [T{c) . The difference 
between Ps{p > Pc) and Ps(pc) is in A and the tail part 
for finite clusters. As p becomes larger than pc, A de- 
creases and the length of tail becomes shorter compared 
to Ps at Pc- The power-law behavior is very robust, be- 
cause it maintains for nearly four decades or 10^ — 10^ 
as shown in Fig. [Ijb) before appearing finite-size effects. 
Moreover, the power-law behavior for p > pc \s nearly 
independent of L as for p < Pc (see Fig. [I{b)). This 
power-law behavior for the finite clusters has been con- 
firmed even for large p upto p = 0.9. In contrast Ps of 
microscopic clusters for p > pc in OP and APCG ex- 
ponentially decays. In 2DSAP the product rule makes 
the macroscopic cluster absorb relatively smaller clusters 
when p gets large in the regime p > pc- Therefore the 
larger microscopic clusters cannot easily disappear. The 
sustainablity of such metastable clusters seems to be the 
origin of the power-law of Pg for p > pc- As we shall see 
the hysteresis of 2DSAP is consistent with the power-law 
for p > Pc because of such metastable states. 

The phase transition for 2DSAP naturally occurs at p 
which divides the two regimes of Pg described by Fig. 
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[Ija) and Fig. [Ijb). The transition threshold pc for 
2DSAP is estimated by the data in Fig. [Ijc) and Fig. 
\V[d). As shown in Fig. [He) and in Fig. [^d) at p = 0.770 
Pg for the macroscopically large cluster starts to be de- 
tached from the continuous distribution of Pg for micro- 
scopic clusters. This detachment behavior seems to be 
independent of L{N) as shown in Fig. [ijc). As shown 
in Fig. md), this detachment behavior barley occurs and 
the hump-like tail still exists for p = 0.765. We have 
scrutinized Pg between 0.765 < p < 0.770, but the sharp 
discrimination between the hump-like behavior and the 
detachment cannot be made. Such complex behavior 
mixing the hump and Ps for the macroscopically large 
clusters for p ^ pc seems to be a unique behavior of 
2DSAP. Therefore the best estimation of Pc from the nu- 
merical data of Pg \spc = 0.768(3). At p ~ pc, Ps satisfies 
the power-law Pg = ^s"* with 5 ~ 0.90(2) very well ex- 
cept for the very tiny tail part. Again this power-law 
Pg = As~'^'^ is very robust and holds for more than four 
decades. The result 5 = 0.9 also provides a very impor- 
tant clue to understand the transition nature of 2DSAP. 
Since Pg is a probability, Pg should satisfy the normal- 
ization condition, Pg = p. However, the summation 
Ps diverges if (5 < 1. Therefore, there should be a 
cutoff Sc in the upper limit as Pi.s) = p. In the 
limit N — > oo, Sc/N — >■ 0. Thus, there should be a dis- 
continuous jump to produce a macroscopic cluster in the 
limit TV — >■ oo and the transition becomes discontinuous. 
The physical origin of the discontinuous transition should 
come from the merge of s-clustcrs with s ~ to form 
the macroscopic cluster when p gets larger to he p = pc- 
One of the most generally accepted and the simplest 
methods to verify whether the observed transition is dis- 
continuous or not is the measurement of the hysteresis 
2l|. The hysteresis measurement for the explosive per- 
colation has also been emphasized in Refs.flll. [l9|. The 
hysteresis is a history-dependent property of a system 
and usually observed in the discontinuous phase transi- 
tion because of the metastable state. If the transition is 
discontinuous, then the route of changes in the order pa- 
rameter Plc during the process of filling sites from the 
vacant lattice or increasing p would be different from that 
for the process of deleting sites from the fully-occupied 
lattice or decreasing p. The order parameter, PlCi is de- 
fined by the probability that a site belongs to the largest 
cluster [l|,i]; 



LC 



Nlc 
N 



(3) 



Here, Nlc is the number of sites in the largest cluster. 
In Fig. [ll^a), to check the existence of hysteresis we com- 
pare the measured Plc's along the process of increasing 
p (solid lines) and along the process of decreasing p from 
N = 2^^ to = 2^°. For the decreasing process, we 
slightly modify the rule (HI) to easily break the larger 
clusters into smaller ones since the rules (II) and 



(a) 1.0 



3- 



o 0.5 



0.0 



- 

L=64 




i=128 




i=256 






L=5\2 




. 1=1024 











0.0 0.2 0.4 0.6 0.8 1.0 

P 



(b) 



0.10 



^ 0.09 
0.08 



200 400 600 8001000 
L 



FIG. 2. (color online) (a) Plot of Plc{p) for the process of 
increasing p (solid line) and that for the process of decreasing 
p (dashed line). L varies from 64 (left most line) to 1024 (right 
most line), (b) Plot of the area, A{L), enclosed by Plc- 



(III) suppress the formation of a large cluster; i.e., if 
tTq > 7r/3 then we delete the site a. With this modified 
rule we find that there exists a hysteresis for various L 
as shown in Fig. HJa). 

Now the remaining question is whether the hysteresis 
robustly remains in the L — > oo limit. For the systematic 
analysis, we show the dependence of area, A{L), enclosed 
by Plc{L) for the increasing and decreasing processes on 
L. If the system undergoes a continuous transition, then 
A{L) should vanish in the limit L — >■ oo. However, our 
data clearly shows that A{L) increases as L increases 
or, at least, seems to saturate to a nonzero value unlike 
the sum rule [l^ in which A{L) — > as L — !> oo. This 
shows that 2DSAP undergoes a discontinuous transition. 
And in 2D lattice, the product rule makes a compl etely 
different transition nature from that of the sum rule |19| . 
This hystcretic property of 2DSAP should be from the 
sustainability of the metastablly larger clusters, which is 
consistent with the analysis of Pg in Fig. [TJ 

Since we don't know the physically corresponding for- 
mula to Hamiltonian or free energy for 2DSAP and there 
exists the nontrivial hysteretic property, it might be phys- 
ically nonsense to discuss about the finite-size scaling. 
However for the comparison's purpose to other works on 
explosive percolation [sl-fiot. we now present the finite 
size analysis around pc = 0.768(3), which is the perco- 
lation transition probability for the p-increasing process. 
From the data in Fig. [DJa), Plc{L) at the Pc is estimated 
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as shown in Fig. [Sja). Plc{L) seems to satisfy the rela- 
tion P(L) ~ L^^ with B = 0.011(2), where convention- 
ally B corresponds to j3/i'. This value of B is very close 
to zero. Thus, in the inset of Fig. ^a), we also fit the 
data to the relation P{L) ~ — logL which corresponds to 
the case i? — )■ 0. Since we cannot exclude the possibility 
_B = or /3 = 0, the possibility for Plc{L oo) at Pc 
to have discontinuous jump cannot be excluded. We also 
measure the mean cluster size, defined by 



S{p,L) = 



(4) 



represents the summation over all s except the largest 
one. S(j),Lys maximal value, Smax{L), is displayed in 
Fig. EJb). Again we fit the data to the conventional 
scaling relation S,nax{L) ~ L^'~^ , and we obtain C — 
1.98(1), where C corresponds to conventional 7/;/. 
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FIG. 3. (a) Plot of Plc{L) measured at pc against L. The 
solid line represents the relation P{L) ~ with B = 0.012. 
Inset: Plot of Plc{L) in semi-log scale. The soUd hne rep- 
resents the relation P{L) ~ — logL. (b) Plot of Smax{L) 
against L. The sohd line denotes Smax ~ L'^ with C — 1.98. 

In summary, we study the site percolation under AP 
with a product rule in a 2D lattice. From the measure- 
ment of Ps{p), we find that the Ps{p) have a very sta- 
ble hump when p < pc. This indicates that below pc 
large number of stable s-clusters with s ~ sjy exist but 
their sizes are still microscopic. As p approaches to Pc, 
Ps {p) has a very robust power-law regime followed by the 
hump. Since the obtained value of the exponent, 5, for 
the power-law regime in the vicinity of Pc is less than 



unity, there should be a cutoff Sc in the possible cluster 
size for p ^ pc unlike OP [l[. Thus, to generate a macro- 
scopic cluster there should be a discontinuous jump in the 
limit L — 00 and the transition becomes discontinuous. 
The non- vanishing hysteresis in P^c also verifies that the 
transition is discontinuous. This result clearly shows that 
the percolation transition caused by the product rule in 
a 2D square lattice is discontinuous. 
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